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Abstract. A generalized WKB approach for constructing WKB series endowed with 
some properties of Taylor ones is presented. Apart from the Riccati equation itself 
its formalism involves also the Riccati-equation's derivatives (REDs) obtained by 
differentiating of the former with respect to a spatial variable. For any smooth potential 
barrier given in the finite spatial interval to include turning points, the zeroth-order 
term of presented WKB series is regular everywhere in this interval. Moreover, the 
more REDs are used, the more exact the zeroth-order solution is. 
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1. Introduction 



In 1994 year two generalized WKB approaches has been proposed - in the paper [Tj, 
for solving the Helmholtz equation, and in the papers [21 [3] (see also [1]), for solving 
the one-dimensional Schrodinger equation (OSE). Both are based on the same idea - 
to use the REDs in order to obtain WKB series regular at the turning points (later, 
in 2002 year, this idea has been reopened in [5]). However, despite the common idea, 
these approaches differ from each other, because they treat the REDs in a different way. 
For example, in the approach [1] the searched-for solutions are expanded traditionally 
in integer powers of Plank's constant h. While in the formalism [2, 3j the role of a 
small parameter is played by some function e(x; h) of the spatial variable x and Plank's 
constant. In this paper we develop the approach [21 H]. 



2. Formalism 



Let us consider the OSE 

^ - Q(°)(z; h)V = 0, QV>(x; K) = ^ (V(x) - E) ; (1) 

where the potential V(x) is zero beyond the interval [a, b] and is a bounded infinitely 
differentiate function; E is the energy of a particle; m e is its mass. It is also assumed 
that there are turning points in the interval [a, b] whose order does not exceed N. 
Remind, if Xq is the turning point of the n-th order (here < n < N), then 

Q(°\x ) = ... = Q^- 1 \x ) = 0; QW(x )^0; (2) 
Q( fc ) = d[Q(o)] K /dx k - if n = 0, points Xq can be formally considered as zeroth-order 
turning points. 

Letting, as in the standard WKB approach, \l/(x) = e s<yX ^ and y(x) = dS{x)/dx 1 we 
reduce solving Eq. (JTJ) to solving the Riccati equation. However, unlike the standard 
WKB approach, we write down this equation together with the first iV REDs: 

y (m+i) + jr cry (m - k) y ik) - Q {m) = 0; ^ + ^ Cf „C"-*>yM - qW = 0; (3) 

k=0 k=0 

m = 0, ... ,N — 1; — dy k /dx k ; C™ are the binomial coefficients. The derivatives 
H) (m = 1, . . . , N) are considered here as independent functions to obey the set 
of Eqs. (j3J); the equality h) = dy( m \x; h)/dx holds, provided that Eqs. ([3]) are 

solved exactly. 
Let now 

yW(x; h) = e~ m ~\x- h)u m (x; h) (m = 0, . . . , N), (4) 

where e(x; h) is some function whose norm \\e(x; h)\\ diminishes when h — ► 0; thereinafter 
the norm ||/(x)|| of any function f(x) is determined as the maximum value of |/| in the 
interval [a, b]. Substitution of Exp. (j3J into Eq. ([3]) yields 

m N 

u m +i + ^2C k n u m ^ k u k -Q m = 0; eu' N - {N + l)e'u N + ^C k u N ^ k u k - Q N = (5) 

fc=0 k=0 



[x;h)\\, \\Q h (x;h)\\ = 0(l), \\u m , k (x;h)\\ = 0(1). (7) 
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where m = 0, ...,N — 1; Qi = e i+2 Q^; £ = 0,...,N; the primes denote the 
differentiation with respect to x. 

Now let us search-for the solution of Eqs. (jSJ) in the form of the expansion in integer 
powers of e(x;h), 

u n (x; h) = ^M m , fc (x; h)e k (x; h), (6) 

fc=0 

supposing that for any k and m the following relations are valid 
d k e(x] H) 
dx k 

Keeping the zeroth-order terms, we obtain 

m N 
ttm+1,0 + ^ C k U m - kfi U kfi = Qm {m = 0, . . . , N - 1) ; Cj?U N - k)0 U kfi = Qn (8) 

fc=0 fc=0 

(we have to stress once more that the functions Qo,- -,Qn are of the same 
order, by norm; it should be taken into account that the function \Qn(x;H)\ takes 
its maximal value at the turning point (of the N-th order), however the ones 
\Qo(x; h)\, . . . , \Qn-i(x; K)\ approach their maxima beyond this point). In the a-th 
order (a = 1, 2, . . .), 

m a 

u m +i, a + ^2 C T u m-k n u Ka ^ = (m = 0, . . . , N - 1); 

k=0 7=0 

du e' N a 

— - (N + 2 - a)-Mjv, a _i ^ u N ^u kta ^ = 0. (9) 

fc=0 7=0 

Now we have to define the unknown function e(x; K) to enter Eqs. ([8]) via the 
functions Q m (x; h) (m = 0, . . . , iV). For this purpose we will suppose that for all m 
u m,o — fn\a m where = 1. Besides, for convenience let further k = e~~ l . Then, from 
Eqs. ([8]), we obtain 

(m) 



k=0 

N 



m + l)a m+ i + V" a m _ k a k = (m = 0, . . . , JV - 1); 

2^ a N- k a k = ^ kN+2 . (10) 

fc=0 

These equations can be reduced to the algebraic equation of the (N + 2)-th order for 
k(x; h). For example, for N = 0, N = 1 and N = 2 we have, respectively, 

k 2 = q(°), k 3 - «g(°) + ^ = o, 3 K 4 -4 K 2 g(°) + «;g( 1 ) + (Q( )) 2 -^ = o. (11) 

3. The relationship between small parameters e(x; ft) and ft at the turning 
point of the n-th order 



Let us now show that there are at least two roots k of Eqs. (ITOl . such that 
limk- 1 ^;^)! =0 (12) 
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for any point x G [a, b\. For example, let x be the turning point of the n-th order; 
\Q^ n+l \x) / Q( n \x)\ <C \k{x\ h)\. The solution of Eqs. ffTUl) . at this point, can be found 
in the form 

oo oo 

k(x;H) = fT^^K s {x)h sP ; a m (x; h) = ^ a m>s (x)H sP (m=l,...,N) (13) 

s=0 s=0 

where (3 is an arbitrary positive constant. The substitution of Exps. ( fT3|) into Eqs. 
( TTTJT) . with taking into account (j2J) and keeping only the main-order terms, yields that 
the series (fT3"j) are justified only for (3 = 2/(n + 2). In this case 

m 

(m + l)a m+ i )0 + ^ a rn-k,oa k ,o = (m = 0, . . . , n - 1, n + 1, . . . , N - 1); 

k=0 

N 

}^ aAr-fc,oQfc,o = (14) 

k=0 

where it is assumed that ao = ao,o = 1- When Eqs. (TT41 have been solved one can find 
Ro from the equation 

2m P 



(n + l)a„ + x,0 + a n-k,0 a k,0 



~n+2 _ ^" b e 



K 



III 



V {n \x ). (15) 



k=0 

From Eq. f|T5|) it follows that there are n + 2 complex roots with equal absolute 
values at the turning point of the n-th order. This means that there are n + 2 different 
functions k(x; h) such that at this point 

\k(x;H)\ ~h- 2/{n+2 \ (16) 

thereby the condition f|T2|) is fulfilled for these roots. If n = 0, i.e., if x is not a turning 
point, then there are two roots to satisfy ( 1T21) . 

So, the expansion in the integer powers of e(x; h) at the turning point of the n-th 
order is equivalent to that in the small parameter h 2 ^ n+2 \ This is in a full agreement 
with the existing approaches. Far from the turning points the expansion is equivalent 
to that in integer powers of h, as in the WKB-approach. 

What is important, one of two relevant roots k(x; H) of Eqs. ffTU]) obeys the condition 
( TT2|) for any point x G [a, b] (hereinafter this root will be denoted as k^'(x; h)). For 
other roots this condition breaks when x crosses turning points. Thus, only one root of 
these equations yields the function e(x; h) to have a small norm in the interval [a, b]. 

The zeroth-order solution, y w [e(x; h)]^ 1 = K^(x;h), associated with this root 
yields a good approximation of the exact solution of the Riccati equation in the whole 
interval [a, b\. In this case, the more the REDs are taken into account, the more precise 
this approximation is. We do not give a strong proof of this statement. Instead we pay 
reader's attention to the fact that for a smooth potential, in the limit iV — ► oo, the 
set of Eqs. ffTUl) for the zeroth order solution coincides with the exact set of Eqs. (j3J), 
considering that in the zeroth order y^ = m\a m K m+l . 

As is seen, the presented WKB series like the Taylor ones includes the derivatives 
of the function expanded; and the larger the number considered derivatives, the more 
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precise the zeroth-order WKB-approximation is. In this sense the presented generalized 
WKB series can be considered as a hybrid Taylor- WKB series. A distinctive feature 
of such series is the possibility to improve the exactness of the zeroth-order WKB- 
approximation. Such possibility is very important, because considering of nonzero terms 
in asymptotic approaches is usually associated with serious mathematical problems (and 
our approach is not an exclusion; however, in our approach, Eqs. (jUJ) are extra in fact). 



4. Example: a linear potential 

To exemplify our approach, let us consider the case when an electron with the energy 
E = O.leV impinges the linear potential V(x) = (l + x/d)Vo to be nonzero in the interval 
[—d/2, d/2]] Vo = O.leV, d = lOOnm. In such setting the point x = is a turning point 
of the first order. 

Figs. Q] and [2] show, respectively, the absolute values of the solutions of the second 
and first equations ffTTj) as well as those of the exact solution of the Riccati equation 
Bi'{px)+i- Ai'(px) 

V{X) = P B i (px) +i .M P x) (17) 
where Ai and Bi are the Airy functions, Ai! and Bi' are their first derivatives, 
respectively; % = \J — 1, p ~ 7,055 • 10~ 2 fira" 1 ; the solution is obtained with the 
help of the Maple programme. Solution ([17]) is evident to correspond to the solution 
*(ar) = Bi(px) + i ■ Ai(px) of Eq. ((H). 

Fig. [H shows that in the limit x — > oo the condition ( fl2|) holds for the roots 1 and 3, 
while in the limit x — > oo it does for the roots 1 and 2 which are the complex conjugates 
to each other in this region. Thus, k^(x; K) is the root 1 for N = 1. For N = 2 (see 
Fig. [2]) the same numeration is used for relevant roots. 

The roots K^(x;h) of Eqs. fflTT) which give the best approximation of y(x) for 
N = 0, iV = 1 and N = 2 can be presented in an analytical form as follows. For 
N = (the first equation in the set (TTTT) ), kW(x;H) = [(Q (0) (z;/i)) 1/2 ]*; hereinafter, 
—ix/n < arg(z l l n ) < n/n for any complex z and integer n; besides, we take into account 
here that if z is a complex root of an algebraic equation with real coefficients, then its 
complex conjugate z* is its root too. 

For N = 1 (the second equation in the set of Eqs. ( TTTT) ) 



h) = X -{A X + A 2 ) - i^(A 1 - A 2 ), (18) 
A 1)2 = surd(r ± \J q 6 + r 2 ,3); q = , r = — ; 

here surd(z, 3) is the Maple's function of a complex variable z: if 9?(z) > then 
surd(z, 3) = z 1 / 3 ; otherwise, surd(z, 3) = — (— z) 1 / 3 ; note that when the function z 1//3 
is used straightforwardly in the Exps. (|18p . this analytical expression yields different 
roots in the different parts of the interval [a, &]. 

From Exps. (TTBl) it follows that far from the turning point x = (when the term 
with QW, in the expression for A12, is negligible) |«;^(a;; ^)| ~ /i -1 ; but at the turning 



Hybrid Taylor- WKB series 



6 




point itself, where Q (0) = 0, = $/\QW/2\ ~ /T 2 / 3 (see also flTSB). 
Lastly, for N = 2 (the third equation in the set (TTTj) ) we have 



«^(a;; ft) = — w(x) + \J (u(x)) 2 — v(x), 



u 




(19) 



«2 



ai 



a 



1 

3" v ~" " v 3 
here z is the solution of the auxiliary equation 



(g(°)(x)) 2 - V), 



X) 



z 3 + b 2 {x)z 2 + h(x)z + bo(x) = 0; 



(20) 



where b 2 = — 02(2), 61 
as follows 

/ \ I / N &2(g) 

2 = si(x) + s 2 (x) — 

6i(x) & 2 (x) 2 



— 4ao(x), 60 = 4ao(a;)a2(x) — (ai(x)) 2 . This solution is defined 



S l,2 



v{x) ± a/ /i(x) 3 + v{x) 2 
b 2 (xf 



1/3 



(21) 



z/ = - [6i(x)6 2 (a;) - 3& (z)] 



3 9 7 6 / \ 27 

Figs. [3] and H] show the ^-dependence of the absolute values and imaginary parts, 
respectively, of k^'(x; h) obtained for iV = 0, A" = 1 and N = 2. As is seen, the more 
the REDs are used in the formalism, the better the zeroth-order approximation is. 

We have to stress that for the potential at hand and for any A" there is always 
such a point x^ in the region x > (the below-barrier region) that S (^'(x; h)) = 
for x > Xi, (see Fig. H|). At the same time for the exact solution y(x) considered here 
there are no regions on the OX-axis where $s(y(x)) = 0. Though |3=(y(a;))| <C |3?(y(x))| 





Figure 3. Absolute values of the exact solution y{x) (circles) and the roots k^(x; H) 
for N = 0, N = 1 and N = 2. 
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Figure 4. Imaginary parts of the exact solution y(x) (circles) and the roots k^\x; K) 
for N = 0, TV = 1 and N = 2. 




Figure 5. Real parts of the exact wave function &(x) (circles) and its zeroth-order 
approximations for N = 0, N = 1 and N = 2. 
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Figure 6. Imaginary parts of the exact wave function ^(x) (circles) and its zeroth- 
order approximations for N = 0, N = 1 and N = 2. 



in the below barrier region, ^s(y(x)) cannot be neglected in any case, because it yields 
the second independent real solution to the OSE, i.e., the function Ai(px). As regards 
the approximate solution k^'(x; h), if Xb £ [a, b] then in the region 6] it loses this 
solution. However, as is seen from Fig. HJ this point shifts to the right on the OX axis 
when N increases. Thus, to eliminate the above shortcoming, one has to include the 
next RED's into the formalism of the zeroth-order approximation. For a given potential 
V(x) and given spatial interval [a, b], there exist such a value of N for which xj, > b. This 
property exemplifies the above statement that the zeroth-order approximation k^'(x; K) 
becomes precise in the limit iV — > oo. 

On Figs. [5] and [6] we show, respectively, the real and imaginary parts of the 
exact solution Bi(px) + iAi(px) of the OSE and its zeroth-order generalized WKB- 
approximations for iV = 0, iV = 1 and N = 2. As is seen, even for the exponentially 
decaying (for x > 0) function Airy Ai(px) the zeroth-order approximation for N = 2 
gives a good exactness. Moreover, the real part of this approximation fits perfectly the 
Airy function Bi(px). As regards the zeroth-order approximation at iV = 0, for this 
case there is no root kP-'{x\ K) which would fit both the independent real solutions of 
the OSE in the below-barrier region {x > 0). 

5. Conclusion 

So, for constructing everywhere regular WKB-series it is suggested to include into the 
WKB formalism the spatial derivatives of a potential-energy function under study. 
WKB-series obtained in such a way are named here the Taylor- WKB ones. For a smooth 
potential given in the finite spatial interval the zeroth-order term of the Taylor- WKB 
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series yields a good approximation in the whole interval under study, including turning 
points if they exist. The more the derivatives of the potential are involved into the 
formalism, the more exact the zeroth-order approximation is. Of importance is the fact 
that for smooth potentials to have only the first- and/or second-order turning points, 
in the considered spatial region, the presented approach yields an analytical expression 
to be the regular approximate complex solution to the one-dimensional Schrodinger 
equation. 
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